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<N . ABSTRACT 

We present the effects of ellipticity of matter distribution in massive halos on 
t"*"- ! the observation of supernovae. A pseudo elliptical Navarro-Frenk- White (NFW) 

mass model is used to calculate the introduced gain factors and observation rates 
of type la supernovae due to the strong lensing. We investigate how and to 
| what extent the ellipticity in mass distribution of the deflecting halos can affect 



surveys looking for cosmologically distant supernovae. We use halo masses of 
1.0 x W 2 h~ l M & and 1.0 x 1O 14 /i" 1 M at redshifts z d = 0.2, z d = 0.5, and 
z d = 1.0, with ellipticities of up to e = 0.2. 
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Subject headings: gravitational lensing — supernovae: general 
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1. Introduction 

X 

Supernovae have emerged as the most promising standard candles. Due to their signifi- 
cant intrinsic brightness and relative ubiquity they can be observed in the local and distant 
universe. Observational efforts to detect high-redshift supernovae have proved their value 
as cosmological probes. The systematic study and observation of these faint supernovae 
(mainly type la) has been utilized to constrain the cosmic expansion history (Goobar & 
Perlmutter 1995; Perlmutter et al. 1999; Schmidt et al. 1998). Light emitted from any ce- 
lestial object is subject to lensing by intervening objects while traversing the large distances 
involved (Kantowski, Vaughan, & Branch 1995) and the farther the light source, the higher 
its chance of being significantly lensed. Apart from the fact that gravitational lensing can 
limit the accuracy of luminosity distance measurements (Perlmutter & Schmidt 2003), it can 
change the observed rate of supernovae as well. 
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Studying supernovae and their rates at high redshifts provide us with much needed 
information for constraining the measurements of the ellusive dark energy, as well as under- 
standing the cosmic star formation rate and metal enrichment at high redshifts. In order 
to observe and, hence, study the faint high-redshift supernovae, one can raise the chance of 
observation by looking through clusters of galaxies or even massive galaxies (see Smail, et al. 
(2002) and the references therein). These 'gravitational telescopes' amplify the high-redshift 
supernovae and thereby increase the chance of their detection. However, this boost in ob- 
servation is offset by the competing effect of depletion (Fig. 1), due to the field being spread 
by the deflector (amplification bias). For an assumed lens model and a given field of view it 
is not obvious which effect dominates the observation of supernovae through the halo. The 
net result depends on the deflector and source parametrs as well as the observational setup 
(Gunnarsson & Goobar 2003). 

Some research has been conducted on the feasibility of observing supernovae through 
cluster of galaxies (see, for instance, Saini, Raychaudhary, & Shchekinov 2002; Gal- Yam, 
Maoz, & Sharon 2002; Gunnarsson & Goobar 2003). These studies have not taken into ac- 
count how the morphology (mainly the ellipticity) of these clusters as gravitational telescopes 
could change the expected supernova rate. In this paper, we investigate whether introduc- 
ing ellipticity into the mass distribution of the deflecting halos can affect the observation of 
supernovae. For this purpose, we use a pseudo elliptical Navarro-Frenk- White (NFW) halo 
model with different values of ellipticity. Throughout the paper we assume the so-called 
concordance cosmology where Q m = 0.3, Q\ = 0.7, and h wo = 0.67, with h wo = H /100 km 
s _1 Mpc _1 . In § 2 we briefly go over the NFW model and show how an analytical formalism 
for a pseudo elliptical NFW mass profile can be introduced. Strong lensing by thin deflectors 
as well as the way ellipticity can afffect the amplification is explained in § 3. We present and 
discuss the results of our calculations in 8 4. 



2. The NFW Halo Model Profile 
2.1. NFW Haloes 

High resolution N-body numerical simulations (Navarro, Frenk, & White 1995, 1996, 
1997) have indicated the existence of a universal density profile for dark matter halos resulting 
from the generic dissipationless collapse of density fluctuations. This density profile does not 
(strongly) depend on the mass of halo, on the power spectrum of initial fluctuations, or on 
the cosmological parameters. These halo models which are formed through hierarchical 
clustering diverge with p oc r _1 near the halo center and behave as p oc r~ 3 in its outer 
regions. Inside the virial radius, this so-called NFW halo profile appears to be a very good 
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description of the mass distribution of objects spanning 9 orders of magnitude in mass: 
ranging from globular clusters to massive galaxy clusters (see Wright & Brainerd (2000) and 
references therein). The NFW halo model is similar to Hernquist profile (Hernquist 1990) 
that gives a good description of elliptical galaxy photometry. However, the two models differ 
significantly at large radii, possibly due to the fact that elliptical galaxies, countrary to the 
dark halos, are relatively isolated systems. 

The spherically symmetric NFW density profile takes the form of 

P(r) = r^\- 2 (1) 
- 1 + - 

r s r s 

where p c = [3H 2 (z)]/(8irG) is the critical density for closure of the universe at the redshift 
z of the halo, H(z) is the Hubble parameter at the same redshift, and G is the universal 
gravity constant. The scale radius r s = r 20 o/c is the charactristic radius of the halo where c 
is a dimensionless number refered to as the concentration parameter, and 

200 

°c = -it— — ( 2 ) 

3 M(l + c)- — 

is a charactristic overdensity for the halo. The virial radius r 20 o is defined as the radius 
inside which the mass density of the halo is equal to 200p c . It is then easy to see that 

M(r 200 ) = M 200 = ^-p c r 3 200 . (3) 

Therefore, NFW halos are defined by two parameters; c, and either r 20 o or M 20 o- For any 
spherical NFW profile with a given mass, the concentration parameter c can be calculated 
using the Fortran 77 code charden.f publicly available on the webpage of Julio Navarro 1 . 

NFW halos can be shown to always produce odd number of images, as opposed to the 
commonly-used singular isothermal sphere (SIS) model which produces either one or two 
images (Schneider, Ehlers, & Falco 1992). Although baryons are expected to isothermalize 
the matter distribution for halos of galaxy mass and below (Kochanek & White 2001), 
taking all of the matter in the universe in isothermal spheres is a great oversimplification 
(Holz 2001). It is, hence, reasonable to model halos (at least massive halos) with NFW mass 
profile instead of SIS model. 



1 http : //pinot.phys.uvic.ca/~jfn/mywebpage/jfnJ.html 
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2.2. Elliptical Potential Model 

Here we present the introduced ellipticity e in the circular lensing potential ip(9), as- 
suming that angular position 9 can be scaled by some scale radius/angle S . The reader 
is encouraged to see Golse & Kneib (2002) and Meneghetti, Bartelmann, & Moscardini 
(2003) for illuminating discussions. We first introduce the dimensionless radial coordinates 
x = (xi,x 2 ) = R/r s = 9/9 s where R is the radial coordinate in the deflector plane, and 
9 S = r s /D(i- Then, one can introduce the ellipticity in the expression of the lens potential 
by substituting x t for x, using the following elliptical coordinate system: 




'Clu Xi 
^0^X2 



= yj xj e + xj e = a u xj + a 2e x 2 2 
= arctan (x2 £ /xi e ) 



(4) 



where a± e and d2 e are the two parameters used to define the ellipticity, as explained below. 



From the elliptical lens potential <p e (x) = (p(x € ), we can calculatete the elliptical deflec- 
tion angle (see § 3.2): 

/ d ^ r \ i— a\ 
= a[x e )sja\ e cos0 e 

(5) 



dxi 



V dx 2 



a(a; e )y / ai7cos e 
Oi{x e ) y/a^ sin e 



/ 



Notice that the expressions above hold for any definition of a u and a 2e - Here, we follow 
Golse & Kneib (2002) who, in order to be able to analytically derive the convergence and 
shear, chose the following elliptical parameters: 



die 
d2e 



1-e 
1 + e 



(6) 
(7) 



which for small values of ellipticity e results in the same ellipticity along the X\ as the 
standard elliptical model of 



ai e = 1 — e 
a 2e = 1/(1 -e) 



(9) 



with e = 1 — b/a, where a and b are the semi-major and semi-minor axis of the projected 
elliptic potential, respectively. 
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3. Gravitational Lensing: a Reminder 



3.1. General Formalism 

In the thin-lens approximation, we define 3 as the optical axis and 3) as the 3- 
dimensional Newtonian potential, with r = \JR 2 + 3 2 . The so-called reduced 2-dimensional 
potential which is defined in the deflector plane is given by 



tpli)) 



2 D 



+00 



ds 



c 2 D d D s 



(10) 



(Schneider, Ehlers, & Falco 1992) where c is the speed of light, and 9 = 9 2 ) is the angular 
position in the image plane. D d , D s , and D ds are angular distances of observer-deflector, 
observer-source, and deflector-source, respectively. The deflection angle a, convergence k 
and the shear 7 are given by the following set of equations: 



f a(9) = V e <p(9) 



k(9)- 



89\ 
i{d)f 



+ 



891 



(11) 



1 / 8 ip 8 ip 



4 \89\ 



891 



+ 



8 2 <p 
89 x 89^ 



The lensing equation then reads: 



(12) 



where (5 = (Pi, P2) is the angular location of the source. The amplification amp of a point 
image formed at 9 is: 

= (1_ K )2_ 7 2 ( 13 ) 

To calculate the angular distances in our work, we use the solution to the Lame equation 
for the distance-redshift equation in a partially filled beam Friedmann-Lemaitre-Robertson- 
Walker (FLRW) cosmology. For a filled-beam flat FLRW cosmology, the angular distance 
D as a function of redshift z is 



D(z) = 



2cz 



1 + z)H (g(z)) 




(^^a) 1/3 ^ 2 
9(z) 



where 



g(z) = 2y/l + tt m z(3 + 3z + z 2 ) + 2 + tt m z(3 + z) . 
See Kantowski (2003) for more detail. 



(14) 



(15) 
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3.2. Lensing Parameters of Spherically symmetric NFW Model 

Several authors have developed the lensing equations for the ordinary, spherical NFW 
halos (e.g. Bartelmann 1996; Wright & Brainerd 2000; Golse & Kneib 2002). Following § 2.2 
we can introduce a dimensionless radial coordinate in the lens plane x = (x 1 , x 2 ) = R/V s = 
9/9 s where 9 S = r s /D d . The surface mass density then becomes 



p(r s x, z)dz = 25 c p c r s F(x) 



(16) 



with 



F(x) = < 



x 2 - 1 
1 



1 , \ _ arcch— 



X 



(x<l) 
(x = l) 



3 

1 ( 1 1\ / 

— 1 . arccos — (x > 1) 

and the mean surface density inside the radius x can be written as 

E(x) 



TlX 



^ / 2irxT;(x)dx = 4<5 c p c r s ^^ 



x. 







with 



i x 1 , 1 / ,\ 

in — | arcch— (x < 1) 

g{x) = l\ + \n- (x = l) 

, X 1 1 / ,\ 

m — | ; arccos — (x > 1 ) 

(see Golse & Kneib (2002)). 

The deflection angle a, convergence k and shear 7 turn out as 



a(x) = 
n(x) = 
7 (x) = 



^crit X 



2«,F(x) 



^crit \ 2< 



= 2k, 



2g(x) 



where k s = <5 c p c r s £ cr - t , with S crit = t 2 D s /(Ai{GD d D ds ). 



(17) 



(19) 



(20) 
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By integrating the deflection angle, the potential (p(x) can be found: 

<p(x) = 2k s 9 2 s h(x) (21) 

with 

In 2 — — arcch 2 — (x<l) 

In — h arccos — (x > 1) 
2 x 



3.3. Lensing Parameters of Pseudo Elliptical NFW Model 

For the particular choice of e in § 2.2, the corresponding convergence and shear can be 
calculated: 



20 s 2 V dx{ dx\ 

e / d 2 ip(x e ) d 2 ip(x e ) 
= K (^) + 7^ 



26 2 V &4 dx 2 t 



K(x e ) + ecos2<f) t '-f(x t ). (23) 



and 



4^ I \ (9x 2 ftr 2 , / \ dxidx 2 ^ 
= 7 2 (x e ) +2ecos20 e 7(f e )/t(x e ) + e 2 (/t 2 (a; e ) -cos 2 20 e 7 2 (x e )). (24) 

Also, the elliptic projected mass density reads: 

S e (x) = S(x e ) + ecos20 e (S(x e ) - £(x e )). (25) 

The lensing equation now becomes (see the appendix): 



Pi = d s x 1 1 - Ak s ei- 



x 2 



= 6 s x 2 [l-4k s e 2 9U 



(26) 



x 2 



and as one expects, the amplification amp reads: 



<*">*>(*) = (1 _ K , (x ) )2 -,|(x) 



(27) 
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It can be shown that ellipticities beyond e = 0.2 result in unrealistic 'peanut' shaped 
projected densities, hence in this work we focus on lower values of e. Figure 2 shows the 
multiple images produced by a 1.0 x 1O 14 /i _1 M halo with ellipticity e = 0.1 (courtesy of 
Golse & Kneib). Dashed lines are the contours with constant surface density S e and the 
solid lines are the critical and caustic lines. Redshifts of source and deflector are 0.2 and 1.0, 
respectively. 

4. The Method 

The main reason for studying supernovae magnified by gravitational lensing is to inves- 
tigate the chance of observing supernovae too faint to be observed in the absence of lensing, 
which is usually the case for cosmologically distant supernovae, specifically type la's. To 
calculate the observed rate of type la supernovae we use the result of predicted rates by 
Dahlen & Fransson (1999) for a hierarchical star formation rate model with a charactristic 
time of r = 1 Gyr (Fig. 3), which limits our calculation to the redshift depth of zuax = 5. 

In order for a supernova to be detected, its apparent magnitude m should not exceed 
the limiting magnitude of the survey mumit- Using the definitions of the apparent magnitude 
and amplification, we get: 

m amp = m + 2.51og(|(l-/€) 2 - 7 2 |) (28) 

in which, vd am ,p is the observed magnitude, and m is the apparent magnitude of the su- 
pernova in the absence of the lensing. We can further write m in terms of the absolute 
magnitude M a t, s of the supernova and rewrite the detection criterion as 

. o on n ( m limit- M abs+ r > \ 

((1 - k) 2 - 7 2 ) D 2 L (z s ) < 10l - ) (29) 

where Dl(z s ) is the luminosity distance of the supernova at redshift z s . The absolute mag- 
nitude of type la SNe has a very narrow Gaussian distribution around M ahs = —19.16 at a 
confidence level of 89% (Richardson et al. 2002) . Here, we assume that the supernova is 
detected as soon as its absolute magnitude becomes brighter than M a b s = —18. 

We take the deflecting halo to be at redshifts z s = 0.2, 0.5, and 1.0, and with virial 
masses of nidi = 1.0 x 10 12 /i~ 1 M Q and m d2 = 1.0 x lO 14 M ft, _1 . Concentration parameter c, 
overdensity 5 C , and virial radius r 2 oo (in units of Kpch^ 1 ) for each case are given in Table 1. 

The field of view is taken to be the spatial angle subtending the virial area of the halo. 
By breaking the projected halo into pixels with the angular size of 5x\ and 5x 2 (which are 
taken to be smaller than the angular resolution of the observation, Figure 4), we calculate 
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the amplification across the halo and hence, find the number of observable supernovae in 
redshift shells with the width of 5z = 0.05. We find the corresponding (spatial angular) 
element S/3i x 5f3 2 in the area behind the halo (in redshift space) where the supernovae are 
bright enough to be detected. Assuming we can arbitrarily minimize 5x\ and 5x 2 , we have 



6(h x 5(3 2 



df3_ 
dx 



Sxi x 5x 2 



(30) 



where 



is the determinant of Jaccobian matrix. 



The reader can refer to the appendix for the derivation of the Jaccobian. The gain 
factor, defined as the ratio of the number of observable lensed supernovae over the number 
of observable supernovae in the absence of lensing (Ni ensed / 'N^oLensing) can be calculated by 
integrating over the predicted rates of both cases across the whole observable area (Fig. 1) 
for any given lensing configuration, considering the ellipticity e. 



5. Results and Discussion 

First, we consider the effect of ellipticity in the number rate of SN la in every redshift 
bin 8z = 0.05. Upper panels of Figures 5 (mai) and 6 (m^) show the number of expected 
supernovae per year occuring in the redshift bins. We present the results for e = 0.0 and 
e = 0.2 with the deflecting halo at redshifts z d = 0.2, 0.5, and 1.0. The survey magnitude is 
assumed to be mu m = 27. The number rate peaks at around z — 1.3 as expected (see Fig. 
3) and dies off rapidly beyond that. It can be seen that the farther the deflector, the slightly 
higher the slope of the curves up to z — 1.3 as a result of higher number of supernovae 
observed in front of the deflector. 

Middle and lower panels in Figures 5 and 6 show the cumulative number rates and the 
gains, respectively. The dominance of amplification bias as a result of the narrowing of the 



m d 


1.0 x 10 12 /i _1 M Q 


1.0 x lO^/i^M© 


Zd 


1*200 


S c 


c 


I"200 


S c 


c 


0.2 


152.61 


38468.6 


9.40 


708.36 


15741.7 


6.46 


0.5 


136.24 


33096.0 


8.83 


632.38 


14426.7 


6.22 


1.0 


111.79 


25118.2 


7.87 


518.88 


12086.5 


5.77 



Table 1: NFW halo parameters for the two halo masses m d at the given redshifts z d used in 
the paper. 
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field in a region immediately behind the deflectors at the assumed redshifts is clear, as the 
gains fall below 1. Beyond that region amplification takes over and more (lensed) supernovae 
are observed. 

In the absence of an intercepting halo, the number rate of the survey drops to zero at 
the redshift limit of the survey. With the deflecting halo present, the observed rate goes to 
zero at a higher redshift. This can be seen in figures 7 (m^i) and 8 (m^) where the deflector 
is at redshift z,i = 0.5 and the survey magnitude limit is mu m = 27. The three upper panels 
depict the expected rates for lensing and no-lensing sccenarios for ellipticities e = 0.0, 0.1, 
and 0.2. The number rates per redshift bin (left) and the cumulative rate (right) are given. 
The reader can readily notice the effect of bias behind the halo. With the galactic size halo 
nidi, the survey can detect supernovae up to redshift z ~ 3 (Fig. 7). This limit increases to 
z ~ 5 (Fig. 8) for the cluster-size halo m d2 . 

The lowest panel in these two figures show the relative difference of the cases with e = 0.1 
and e = 0.2 with respect to e = 0.0. The 2 curves do not show significant difference for the 
redshift bins in front of the halo. In the regime behind the halo, the difference becomes 
remarkable: it increases up to redshift z = 1.4 for m dl and z = 1.7 for m d2 . The difference 
doesn't vary remarkably beyond the maximum point. 

To further see how ellipticity changes the expected rate of observed supernovae we put 
the result of our calculations for different ellipticities for a given range of magnitude limits 
on the same plot. Figure 9 shows the number rate of observed type la supernovae (upper 
panel) for ellipticities e = 0.1 and e = 0.2 together with their relative difference with respect 
to the case with no ellipticity (lower panel). Both halo masses, m d i and m d 2 are at redshift 
z d = 0.2. Figures 10 and 11 show the results of the same calculations with halos at redshifts 
z = 0.5 and z = 1.0, respectively. The number rates in each figure increase smoothly up to 
the magnitude limit at which the survey is deep enough to detect the supernovae as far as the 
halo itself, e.g, mu m = 22.4 for a concordance cosmology of (fl m ,u\, foioo) = (0.3, 0.7, 0.67). 
From that point on the rates increase very rapidly as the magnitude limit goes up. That is 
caused by the halo lensing and hence amplifying the supernovae which would otherwise be 
too dim to be observed. The relative differences depicted in these figures show that even 
at a magnitude limit of 25, effect of ellipticity cannot be ignored as it significantly changes 
the number /percentage of the observed supernovae; for instance, the relative difference for 
e = 0.2 with deflecting halo m^ 2 at redshift z = 1.0 (Fig. 11) exceeds 9% for the magnitude 
limit of mn m = 27. 
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6. Conclusion 

Aiming behind massive halos seem to be a good way to enhance the high-redshift su- 
pernovae surveys. The cumulative gains of such surveys seem insignificant at low redshifts 
(z s < 0.2) but the results are remarkable at higher redshifts. For deep observations where 
mn m > 25, the geometry of the intervening halo cannot be ignored. We have shown that 
introducing ellipticity in the (gravitational potential of) the mass distribution of a deflecting 
halo (here, for a galactic halo of mass 1.0 x 10 12 M Q /i _1 as well as a middle-size cluster of 
galaxies with a mass of 1.0 x 10 14 M Q /i -1 ) can affect the rate of observed supernovae by a 
few percent. It was shown that the farther the supernova survey probes, the more significant 
the effects of introduced ellipticity become. 

It should be noted that this work does not involve a broad range of mass profiles for the 
halos (although we specify that the survey is limited to the virial area of the halo), nor does 
it address the much needed k-correction. Our calculations are actually an oversimplification 
due to the fact that a large, massive halo like a galaxy cluster has substructure which consists 
of the member galaxies, as well as large clouds of gas. A more sophisticated lens model with 
ellipticity should be employed to calculate the number rate of observed supernovae. 

The author wishes to thank D. Branch for enlightening discussions and suggestions. 
The author would also like to thank C. Gunnarsson and A. Goobar for generously offer- 
ing him their data set on the rate of SNe la, and J. Navarro for allowing him to use the 
NFW code. This work was in part supported by NSF grant AST0204771 and NASA grant 
NNG04GD36G. 



A. Appendix 

Here we derive the lensing equation for an elliptical NFW halo with ellipticity of e 
introduced in its 2-dimensional potential, and proceed to calculate Jaccobian ^ needed to 
get the spatial angular element 5f3\ x 5f3 2 in the source frame. 

Lensing Equation 

Introducing the dimensionless coordinate system x = (xi,x 2 ) = R/r, 
lensing equation becomes 

f A = °sXl - Ct\ [Xl,X 2 ] 

\ (3 2 = s x 2 - a 2 [x 1 ,x 2 ] 



= 0/0, , the 
(Al) 
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Given the elliptical deflection angle of 

a e (x) = 



a(£ e )v/ai7cos0 ( 
y -q— = a{x e )y/a 2e smtp € 



and the deflection angle of a as 



a(x) = 6— — = 4K s —g(x)e x 

^crit X 



the lensing equation now reads 



#2 = 6>s^2 1 - 4A; s e 2 



(A2) 



(A3) 



(A4) 



Jaccobian 



To calculate spatial angular element 5(3i x <5/3 2 we use the Jaccobian equation 



5f3 l x 5(3 2 
Given Eq. A4, we get: 

r d(h 

dxi 
dxi 
dx 2 
dx 2 



dp 



dx 



5xi x 5x 2 



dp, dp 2 dp, d(3 2 



dx, dx 2 dx 2 dx. 



8x1 x 5x 2 



6 S (1 - 4k s a u G(x e )) 9 sXl (l - 4k s a le ^^ 
dG(x e )' 



6 s x 2 1 - Ak s a 2t 



9 s xi ( 1 — 4k s a, ( 



dx, 
dG(x £ )' 

dxo 



= 6 s {l-Ak s a 2e G{x e ))6 s x 2 [l-Ak s a 2 



dG{x e 
dx 2 



where function G is defined as 



and 



x f 



(1 -x?)' 
(1 + ^ 



xi 



arcch 



xi 



x e 2x e (1 — x\ 



x f 



2x e (l-x 2 e ) (l-x 2 e ) 



arccos- 



x f 



(X t < 1) 

(X e = 1) 
(X t > 1). 



(A5) 



(A6) 



(A7) 



(A8) 
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Fig. I.- 

Schematic picture of the lensing configuration by a deflecting halo. Zhaio is the redshift 
of the halo and zi imit is the redshift corresponding to the limiting magnitude m^t. The 
shaded area shows the volume where SNe are bright enough to be observed. 

Fig. 2.- 

Multiple images produced by a 1.0 x lO 14 /^ 1 M halo with ellipticity e = 0.1. Dashed 
lines are the contours with constant surface density and the solid lines are the critical and 
caustic lines. Redshifts of source and deflector are 0.2 and 1.0, respectively (courtesy of 
Golse & Kneib). 

Fig. 3.- 

Rates of type la supernovae in intervals of Sz = 0.05. Dilution factor of 1 + z is taken 
into account (courtesy of Gunnarsson & Goobar). 

Fig. 4.— 

This figure shows how the projected deflector is 'pixellated' in order to calculate the 
observable area behind the halo. Each pixel has dimensions of 5u x 5uj with 5uj being smaller 
than the angular resolution of the observation. 

Fig. 5.- 

Rates of observed supernovae la per redshift bin S z = 0.05 (upper panel), cumulative 
rate (middle panel), and the lensing gain (lower panel) for a deflecting halo of mass rrid = 
1.0 x lO^/i^M© at redshifts of z d = 0.2, 0.5, and 1.0 with ellipticities e = 0.1 and 0.2. 

Fig. 6.- 

Same as Figure 5, with m d = 1.0 x 10 14 M /i -1 . 
Fig. 7.- 

In this picture the 3 upper panels show observed rates of lensed (solid line) and unlensed 
(dash line) for three different ellipticies e = 0.0, 0.1, and 0.2. The deflecting halo has a mass 
of rrid = 1.0 x 10 12 /i -1 M Q and is located at redshift z s = 0.5. The lowermost panel depicts 
the relative difference of e = 0.1 (solid line) and e = 0.2 (dash line) with respect to e = 0.0. 



Fig. 8.- 

Same as Figure 7, with m d = 1.0 x lO 14 M /i _1 . 
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Fig. 9- 

Rates of observed supernovae la as a function of survey magnitude limit m (upper panel). 
Result sare shown for halo masses m d = 1.0 x 10 12 /i~ 1 M Q and m d = 1.0 x lO u h~ 1 M Q with 
ellipticities e = 0.1 and 0.2. The halo is at redshift z d = 0.2. The relative difference of cases 
with e = 0.1 and 0.2 with respect to e = 0.0 is given in the lower panel. 

Fig. 10 — 

Same as Fig. 9 with z d = 0.5. 



Fig. 11.— 

Same as Fig. 9 with z d = 1.0. 
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